arXiv:hep-ph/9808477vl 30 Aug 1998 


Sine-Gordon Parametric Resonance 


Patrick B. Greene 

Department of Physics, University of Toronto, 60 St George Str, Toronto, ON M5S 1A7, Canada * 

Lev Kofman 

CITA, University of Toronto, 60 St George Str, Toronto, ON M5S 1A1, Canada I 

Alexei A. Starobinsky 

Landau Institute for Theoretical Physics, Kosygina St. 2, Moscow 11733f, Russia 

(February 1, 2008) 

We consider the instability of fluctuations in an oscillating scalar field which obeys the Sine-Gordon 
equation. We present simple closed-form analytic solutions describing the parametric resonance 
in the Sine-Gordon model. The structure of the resonance differs from that obtained with the 
Mathieu equation which is usually derived with the small angle approximation to the equation for 
fluctuations. The results are applied to axion cosmology, where the oscillations of the classical axion 
field, with a Sine-Gordon self-interaction potential, constitute the cold dark matter of the universe. 
When the axion misalignment angle at the QCD epoch, 9 o, is small, the parametric resonance of 
the axion fluctuations is not significant. However, in regions of larger 9q where axion miniclusters 
would form, the resonance may be important. As a result, axion miniclusters may disintegrate into 
finer, denser clumps. We also apply the theory of Sine-Gordon parametric resonance to reheating 
in the Natural Inflation scenario. The decay of the inflaton field due to the self-interaction alone 
is ineffective, but a coupling to other bosons can lead to preheating in the broad resonance regime. 
Together with the preheating of fermions, this can alter the reheating scenario for Natural Inflation. 
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I. INTRODUCTION 

An oscillating pseudo-Nambu-Goldstone (pNG) scalar 
field, (j>, governed by the Sine-Gordon theory and the 
parametric resonance of its fluctuations, 4>ke lkx , is a pro¬ 
totype for several interesting cosmological and other ap¬ 
plications. One of these is the cosmic axion. The axion 
field, which was originally invented in order to solve the 
strong CP problem ||], is known to be a dark matter can¬ 
didate @1 , for reviews see e.g. 10 - In the literature, it 
is often assumed that the axion field is a relatively homo¬ 
geneous, oscillating classical scalar field which does not 
decay to other particles. The possibility of a paramet¬ 
ric resonance decay of this field due to its self-interaction 
was considered in Ref. Using the small angle ap¬ 
proximation, which leads to a Mathieu equation, and 
an elementary estimate for the width of the instability 
band and the redshifting of the modes due to the expan¬ 
sion of the universe, it was concluded that this process 
is inefficient. Another important application of the Sine- 
Gordon theory is to Natural Inflation Q, where after 
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inflation a pseudo-Goldstone field oscillates coherently 
about the minimum of its Sine-Gordon potential. The 
parametric resonance of inflaton fluctuations due to the 
self-interaction in the natural inflation scenario was con¬ 
sidered in ||. Again, using the small angle approxima¬ 
tion and the Mathieu equation, it was concluded that 
the resonance is not effective. One more application is 
the resonant decay of disordered chiral condensates [^o| , 
modeled with the Sine-Gordon equation [^lj . 

In the last few years, it has been appreciated that para¬ 
metric resonance of coherently oscillating scalar fields, 
in particular, the inflaton field (preheating) can be very 
efficient [|l2| . For this reason, axion parametric res¬ 
onance was recently re-examined ]l^ |. The equations 
governing axion fluctuations were solved numerically in 
Minkowski space-time and in an expanding universe un¬ 
der the simplifying assumption that the axion mass is a 
non-vanishing constant since the very beginning of the 
radiation dominated universe. The results confirmed the 
conventional wisdom that axion parametric resonance is 
inefficient because of the expansion of the universe. How¬ 
ever, in the cosmic axion scenario, the axion mass is 
initially absent; it turns on very abruptly at the QCD 
epoch, tQCD , when instanton effects break the classical 
Peccei-Quinn symmetry, Upq( 1). After Iqcd , the axion 
begins to oscillate with a period many orders of magni¬ 
tude smaller than the age of the universe. This is differ¬ 
ent from the setting of ref. 1131, where consequently, the 
effect of expansion is overestimated. 
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In this paper we will systematically study the para¬ 
metric resonance in the Sine-Gordon theory, both in 
Minkowski space-time and in an expanding universe. We 
will show how this problem can be solved analytically 
and present the theory of the parametric resonance of 
the Sine-Gordon equation. For this, we will use meth¬ 
ods which we elaborated in application to the theory of 
preheating after inflation (ljjl^] . With these novel tools, 
we will reexamine the axion parametric resonance in an 
expanding universe where the axion mass is switched on 
at tQCD- It turns out that the results of the Sine-Gordon 
parametric resonance can be extended to axions in an ex¬ 
panding universe analytically. In agreement with j~i||7Ti| 
we will show that the axion parametric resonance is in¬ 
effective as soon as the axion misalignment angle at the 
QCD epoch, 8q, is small, 0 o ~lO -3 . 

However, there are axion scenarios where the axion 
field 8q is spatially inhomogeneous, and in its peak re¬ 
gions can be significantly larger than its mean ampli¬ 
tude. When cosmic axions become gravitationally domi¬ 
nant, these regions will turn into dense axion miniclusters 
fiT|pT| . We consider the axion resonance in regions where 
the coherently oscillating axion field of the large ampli¬ 
tude 9q is smooth at the scales of proto-miniclusters. 
In these regions, the axion parametric resonance can be 
quite significant and leads to the decay of the smooth 
axion field into an inhomogeneous field within proto- 
miniclusters. Then, when axions become gravitation¬ 
ally dominant, these inhomogeneities will form small high 
density clumps inside the minicluster. 

Finally, we turn to another application of the Sine- 
Gordon parametric resonance: reheating in natural in¬ 
flation. We consider the parametric resonance of inflaton 
fluctuations due to the self-interaction. We confirm the 
result of although by different reasoning, that para¬ 
metric resonance due to the self-interaction is not effec¬ 
tive. For the sake of completeness, we also consider the 
creation of other bosons x via a coupling of the form 
g 2 (f> 2 X 2 with the inflaton. In this case, X"P ar ticles can be 
created in the broad parametric resonance regime; and 
thus, preheating in natural inflation can be very efficient. 


II. PARAMETRIC RESONANCE IN THE 
SINE-GORDON THEORY 

First, we neglect expansion of the universe (a = 1, 
H = ^ = 0), and consider the parametric resonance in 
the Sine-Gordon theory. We will start with the time evo¬ 
lution of the homogeneous background oscillations in this 
theory. Then we turn to the basic issue of our study, the 
stability of small inhomogeneous fluctuations 
of comoving momentum k. Since the background oscilla¬ 
tions are periodic, the temporal part of the eigenfunction 
is 4>k(t) ~ e /ifc7rat , where the characteristic exponent ^ is 
real for resonant modes. The goal of this Section is to 
find fik in the Sine-Gordon theory. Fortunately, one can 


avoid making any further simplification (besides neglect¬ 
ing expansion) because it is possible to find exact reso¬ 
nance solutions using the methods developed in Ref. [^5| . 
Once we have developed the solutions, we compare the 
exact result with the estimation obtained from the small 
angle approximation, 9q <C 1, which leads to a Mathieu 
equation for the fluctuations. 


A. Background Oscillations in the Sine-Gordon 
Theory 

The potential of the scalar field in the Sine-Gordon 
theory can be represented as follows: 

V(<j>) = A 4 (l-cos , (1) 

where A 4 = m 2 f 2 with / a constant (corresponding to 
the radius of the tilted “wine bottle” potential) and m the 
mass of the bose field q !>, see Fig. [j]. To simplify notation, 
we will sometimes use the dimensionless misalignment 
angle, 9 = j. 



± 

f 

FIG. 1. The potential A 4 ^1 — cos j'j . 

In general, the field </>(f, x) obeys the Sine-Gordon 
equation. For a homogeneous scalar field cf> = </>(f) and 
without expansion of the universe, we have 

(j) + m 2 f siny = 0 , (2) 

where m is a constant mass and / is the radius of the po¬ 
tential. For small initial misalignment angle, 9q = y- -C 
1, the solutions to this equation are approximately har¬ 
monic. Nevertheless, we will derive the general solution 
of this equation for an arbitrary initial amplitude without 
the assumption of a small misalignment angle. Although, 
for practical purposes in dealing with axion parametric 
resonance one may skip this general consideration, it is 
instructive to consider parametric resonance in the Sine- 
Gordon theory for a wide range of parameters. This will 
eventually allow us to understand how the parametric 
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resonance in an expanding universe can be so sensitive to 
the initial conditions and parameters involved. 

Without loss of generality, we can restrict ourselves to 
the finite motion \<j>(t)\ < 7r/. The energy integral of 
Eq. (g) is 

t + V ( ( / ) ) = E, (3) 

where E < 2m 2 / 2 , V((j>) is given by Eq. (g), and E = 

v(M. 

In Appendix A we derive the final result for the evolu¬ 
tion of the background oscillations 


(j>(t) = 2/ arctan — cn(mt ; /C)J , (4) 

where the useful dimensionless parameter is e = = 

1 — cos Oq and the modulus K, = y = \f\- As follows 
from Eq. (g), it is convenient to use the dimensionless 
time variable 

r = mt. (5) 

As shown in Appendix A, the period of the background 
fluctuations in the dimensionless time variable, T t , can 
be expressed through a complete elliptic integral of the 
first kind: T t = 4K^y'A). One may check that the 

amplitude of <j>(t) from Eq. (g) is <j> o- 

In the small amplitude limit, (j >0 <C /, we have E sa 
2 , 2 

J >0 and the oscillations of the background field are ap¬ 
proximately harmonic: <p(t) ~ c/>o cos mt, T t « 27t. 

B. Analytic Resonant Solutions 

We shall consider the fluctuations (f> k (t) in the Sine- 
Gordon model without expansion of the universe, but 
not assuming the misalignment angle -y- to be small. In 
this case the mode equation for fluctuations is 

4>k + ^ k 2 + m 2 cos j ^4> k = 0 . (6) 

Let us use the dimensionless time, r = mt, and the di¬ 
mensionless momentum k = —. Then the mode equation 

m 1 

is 

4>k + p 2 + cos <t>k = 0 , (7) 

where the time-dependence of the background field, </(r), 
is given by Eq. Recall that </>(r) is a periodic function 
of t with period T t . 

Equation (0) describes an oscillator, ^(r), with a vari¬ 
able frequency which depends periodically on time. As 
is well known, the solutions 4>k{r) can be exponentially 


unstable, (j>k oc e ^ kT , for some ranges of the parameters 
k and y . For positive-frequency vacuum initial condi¬ 
tions, ~ y=y, one expects the exponentially fast 

creation of particles nk oc e 2MfcT in the unstable modes 
as the background field oscillates. For a given mode, k, 
the strength of interaction is determined by the dimen¬ 
sionless background amplitude, or misalignment angle, 
do = y. That is, the absolute value of the misalignment 
angle, 9q, ultimately defines the structure of the paramet¬ 
ric resonance in Sine-Gordon theory. It turns out that the 
strength of the resonance monotonically depends on this 
parameter; it is stronger for larger |0 q|- For |0 q| > 7t/ 2 
the effective mass is negative, and one has even a strong 
tachyonic instability. 

The main purpose of our study in this Section is to find 
the characteristic exponent pk = Pkifi o)- Even neglecting 
the expansion of the universe, the equation for fluctua¬ 
tions (g) looks difficult at first glance. Fortunately, it 
turns out that in the Sine-Gordon theory one can ob¬ 
tain simple, closed-form analytic solutions to the mode 
equation (g). 

For this, let us rewrite this equation in a different form. 
We will use a new “time” variable z : 

z(t) = cos P = —\/2(l - z 2 )(e - 1 + z)~P , 

( 8 ) 

where, as before, e = ^ = 1 — cos^o- Note that z is 
sandwiched in the range 1 — e < z < 1. Equation (0) for 
fluctuations becomes 

2(1 — z 2 )(e — 1 + z)<j>k + ^1 + 2(1 — e)z — 3z 2 ^0' k 

+(k 2 + z)(f>k = 0 , (9) 

where (...)' stands for the ^-derivative. As we show in 
Appendix B, the mode equation in this form coincides 
with the algebraic form of the Lame equation with par¬ 
ticular coefficients. 

At this point we will repeat the trick we found in 
for a similar problem in A</ 4 theory. Specifically, we intro¬ 
duce two linearly-independent solutions of Eq. ®), Mz) 
and (j >2 {z), where the lower index k is omitted for sim¬ 
plicity. Let us construct the bilinear combinations </> 2 , 
<j> 2 , and 4>i4>2- From (If) it follows that these bilinear 
combinations obey a third order equation 

2(z 2 - l)(e - 1 + z)M"' + [9z 2 - 6(1 - e)z - 3] M" 

+2 {z - 1 + e - 2k 2 )M' - 2M = 0 . (10) 

The three solutions, M(z), of this equation correspond 
to (up to a numerical normalization A/ still to be given) 
the three bilinear combinations of (f> 1 and 02 - 

Equation ( |l0|) admits a polynomial solution 

M{z) = z — 1 + e — 2 k 2 . (11) 
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In the resonance zone, this polynomial solution must be 
the product of an exponentially growing solution and an 
exponentially decreasing one, i.e. 

4>i(z)4>2{z) = J\T 2 ■ (z - 1 + e - 2k 2 ) , (12) 


where J\f is a normalization factor. An additional relation 
between the two fundamental solutions is given by the 
Wronskian of Eq. (||) 


01^2 - ^2^1 = A/ -2 


C-k 


V(l-* 2 )(e-l + *) ’ 


(13) 


where Ck is a constant to be determined. From © and 
© we obtain the closed form analytic solutions 


<t>iA z ) = A/V|M(z)| x 

-iA^wi ) ^ (I4> 

The solution to Eq. ( 0 ) with positive-frequency vacuum 
initial condition will be a linear superposition of the two 
fundamental solutions <f>i .2 jn|. The amplitude of fluc¬ 
tuations is rapidly dominated by the growing term. As 
is shown in Appendix C, we can extract the character¬ 
istic exponent describing the resonance effect from the 
fundamental solutions themselves. Therefore we will not 
further discuss the numerical normalization and, for sim¬ 
plicity, we normalize our fundamental solutions to unit 
initial amplitude: A f = |M(1)| -1 / 2 . 


C. The Width and Strength of the Resonance 


The closed-form analytic resonant solution ([hi]) allows 
us to find the structure of the resonance in terms of the 
width of the resonance band Are 2 and its strength Hk for 
an arbitrary choice of parameters. Indeed, substituting 
this solution back into equation (||) for <j>(z), we find the 
constant c/.: 

cl = 8 k 2 (| - k 2 ) (f - | + k 2 ) • (15) 


For exponentially growing solutions, Ck must be real. 
Therefore, the exponentially growing solutions for fluc¬ 
tuations with re 2 > 0 occur in a single instability band 
for which 0 < re 2 < |. In terms of the basic parameters, 
the width of the instability band is 


0 < re 2 < 


1 — cos d 0 
2 


(16) 


The growing solution of Eq. ( 0 ) has the form <pk{ T ) = 
e^ kT P[ z{t)\, where P\z(t)\ is a periodic function of time 
r. Using Eq. (p~4|), we can now find the characteristic 
exponent /r*, as a function of re for a given initial mis¬ 
alignment angle, do- The technical details can be found 
in Appendix C. The result is 


Mfc(0 o ) = \j 2re 2 ^l — (cos d 0 + 2re 2 ) 2 ^ , (17) 

where an auxiliary function J K = I k (0q) is 

7r/2 

4(00) =4 

sin2 19 _ (18) 
1 + (cos do + 2re 2 ) sin 2 id 

Recall that T T (do) = 4K^|sin^p|^ is the period of the 
background oscillations. 

In summary, we have found that for an arbitrary value 
of 9q there is always a single resonance band given by for¬ 
mula &■ The characteristic exponent strongly depends 
on the amplitude do and its value is given by an explicit 
integral in formula ©)■ The characteristic exponent /ifc 
found by numerical integration of the equation for fluc¬ 
tuations © with vacuum positive-frequency initial con¬ 
ditions is given in Fig. | as a function of re for various 
values of the initial background misalignment angle do- 
These results coincide precisely with the analytic formula 

©• 


did 


yj (1 + sin 2 $)(1 + cos do sin 2 id) 



FIG. 2. The characteristic exponent /r*, of the 

growth of fluctuations (<j>k oc e MfcT ) as a function of 
the (scaled) resonant momentum re in Sine-Gordon the¬ 
ory. Different curves correspond to different values of 
the initial background misalignment angle, do. From 
the outer curve inwards do = 1.67,1.25,1.0,0.83,0.71, 

0.625, 0.56, 0.5, 0.45, 0.42, 0.38, 0.36. 


Let us now consider the case of small misalignment 
angle do <C 1, which is very important for applications. 
In the small angle limit, the general formula (hd) is re¬ 
duced to the very simple expression fi.k = f \J$o ~ 4k 2 ; 
or explicitly in terms of the momentum fc, 



(19) 


Q 2 ,2 

In this case, the maximum value Uk = occurs at -A = 

’ 16 m z 

Q 2 ,2 $2 

The resonance band in this limit is 0 < Ay < In 
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the next Section we will extend formula (|lt]) taking into 
account the expansion of the universe. In this form it 
will be used to study axion parametric resonance, where 
the background axion misalignment angle will be small 
by the QCD time, tQCD- 


D. On the Validity of Use of the Mathieu Equation 


Let us return to the equation for fluctuations (0) and 
consider a small misalignment angle, 9 = j <C 1. There 
we will make the approximation: cos 9(r) ss 1 — ^ 9 2 (t ) ss 
1 — ^9 2 cos 2 r. Thus, to this order in 9 2 , Eq. © becomes 


</>fc 


+ 



9 2 

-j- cos (2r) 
4 


4>k = 0 , 


( 20 ) 


which is a form of the Mathieu equation, cj>k + 
{A — 2gcos 2r) </>*, = 0, with parameters A = 1 — 2q + k 2 

q2 

and (/ = j < 1. The properties of the Mathieu equation 
for small q parameter are well known 1171. The resonance 
in the first (leading) band occurs for A = 1 ± q, which 

q 2 30 ^ 

gives a width for the resonance band -§- < k 2 < -g 2 -. 
The characteristic exponent in the first zone is fj,jf = 

\j (f) 2 — (V^A — l) 2 or, explicitly to this order, 


»k = 




( 21 ) 


Note that this expression differs from estimations in 
H. In the Mathieu approximation the leading resonance 

band is concentrated at k 2 ~ -j-m 2 . 

We presented here this analysis because it is simple 
and intuitive. The approximate solution obtained with 
the Mathieu equation has some qualitative features simi¬ 
lar to the exact solution, but it is quantitatively different. 
In particular, the actual resonance band begins at k = 0. 
Even though this band may be very narrow, the red- 
shifted resonant modes <pk never shift below this band, 
i.e. expansion of the universe is not as destructive as 
one might expect. In this respect, the actual resonance 
solution is qualitatively different from the approximation 
based on the Mathieu equation. 


III. BACKGROUND AXION FIELD 

In this and the next section we will investigate ax¬ 
ion parametric resonance using the results derived above. 
Before considering axion fluctuations, we shall discuss the 
evolution of the background axion field. In the context 
of the problem of fluctuations, it can be either homo¬ 
geneous or quasi-homogeneous depending on the cosmic 
axion scenario and relation of the length scales of fluctu¬ 
ations to the scale where the background field is smooth. 


The axion, q. i, is a pseudo-Nambu-Goldstone boson 
which appears after spontaneous symmetry breaking in 
a theory of a complex scalar field. It acquires a small 
mass, m a , due to instanton effects when the temperature 
of the universe drops below the QCD scale @0- The ef¬ 
fective potential in the simplest theory of the axion field 
can be represented as ([j]), where / is the radius of the 
axion potential [lj. The axion is massless at very high 
temperature but gradually develops a mass. For T A 
we have m(T) « 0.1m a (A/T) 3 - 7 . Around the QCD phase 
transition (T ~ A ss 200 MeV), the mass is quickly sat¬ 
urated at its low-temperature value, m a = (A 2 //), or 

m a « 0.6 x 10“ 5 eV. 

At this point, the further discussion of the background 
axion field branches into different axion scenarios: the 
initial axion misalignment is globally homogeneous, the 
initial misalignment is homogeneous within a Hubble do¬ 
main at the ‘defrost’ moment 1 1 , or the axions are pro¬ 
duced by the decay of axionic strings and walls. To in¬ 
vestigate the axion parametric resonance due to the self¬ 
interaction of the fluctuations with the background axion 
oscillations, in principle, one needs to specify the back¬ 
ground axion field <p(t, x) in these models. However, as we 
will see below, the axion parametric resonance depends 
on the value of 9q at the moment of the QCD phase tran¬ 
sition. As soon as a region with given misalignment angle 
9q is smooth and much larger in size than the typical res¬ 
onant wavelength, one can split the axion field <p(t, x) 
into a background and fluctuations around it. 

We first consider the simple model where cf>(t) is a co¬ 
herent homogeneous field throughout the universe, and 
then will discuss the inhomogeneous background field 
<j>(t,x). A homogeneous axion field obeys the Sine- 
Gordon equation in an expanding universe 

4> + 3H^ + m 2 {T)f sin j = 0 . (22) 

The time evolution of the background axion field is plot¬ 
ted in Fig. |j. Before the condition m(T) « 3 H(t) holds 
at a moment t\, the axion field is frozen at its initial value 
<f >i (corresponding to a misalignment angle 9\). Shortly 
after this ‘defrost’ moment t±, the field begins to oscil¬ 
late with gradually decreasing amplitude. After the axion 
mass is saturated at tQco-, the coherent fluctuations are 
approximately harmonic, ~ </>o > where </>o is 

the amplitude at tQCD , and a = \JtjtQCD■ 
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FIG. 3. Time evolution of the background field 0(t), in 
terms of the misalignment angle 9(t) = — ^, obtained with 
the numerical solution of Eq. (^). The initial value is 
chosen, for illustration, as the r.m.s. of the random phase, 

i.e. do = 

The resonance regime for axion fluctuations starts at 
tQCD (the growth of fluctuations is negligible prior to 
that moment). Therefore, we have to evaluate the mis¬ 
alignment angle 9q at tQCD- If the mass density of co¬ 
herently oscillating axions in units of the critical density 
is fl a , then 

el ~ ict 6 ^ 2 , (23) 

where h is the present Hubble parameter in units of 
100 km/s • Mpc. We assume a temperature at 
Tqcd — A qcd — 200 MeV. The value given by Eq. ( j23| ) 
is very small and this, we will see, suppresses the axion 
resonance. 

There are, however, cosmic axion scenarios in which 
the axion field at the QCD epoch is globally inhomoge¬ 
neous. In this case, the free axion field (f>(t , x) obeys the 
Sine-Gordon equation (|2^) with an additional gradient 
term 4j-V| (j>. In the uncorrelated misalignment scenario 
fj[[i there are rare (overdensity) domains where the ini¬ 
tial phase 6*i is as large as n. As was shown in [p~8| , 
large inhomogeneities of the background axion field at 
the present-day scales t\ X (1 + z\) ~ 10 18 cm upwards 
result in large density fluctuations by the epoch when 
axions gravitationally dominate. This leads to the for¬ 
mation of axion miniclusters of the mass ~ 1O _12 M 0 . 
For our purpose, it will be important that the misalign¬ 
ment angle at tQCD i n the overdensity domains can be 
significantly larger than the average given by (f23[). The 
size, A, of these domains, or proto-miniclusters , at tQCD 
is given by the size of the horizon at the defrost moment, 
fi, stretched by the scale factor a(tQCD)/a(t i). This 
size L is smaller than the size of the horizon at tQCD , 
but significantly larger that the wavelengths ~ (#o?7i a ) -1 
of axion fluctuations subject to resonant amplification. 
Therefore, one can consider the resonance of axion fluc¬ 
tuations around a (locally) homogeneous and coherently 
oscillating field 9 with an amplitude considerably greater 
than the estimate (|23|). 


The axion field is also inhomogeneous in the scenario 
where axions are radiated by axionic cosmic strings @ 0 - 
However, the history of 9(t) depicted in Fig. is not 
applicable to this scenario. Again, the mean value of 9q 
is given by formula (|2^), but the axion field at tQCD is 
inhomogeneous. Let us briefly discuss the properties of 
9o(x) in this scenario. 

In this scenario, according to Refs. [Q, the axion field 
appears after the universe cools down and PQ symme¬ 
try breaking occurs. In different parts of the universe, 
on the order of the size of the horizon at that moment, 
the axion field takes different values in the range from 
0 to 2jrf. As a result, an axionic string network forms. 

Starting with the moment f* ~ 10~ 12 sec, the 

strings begin to oscillate and radiate their energy into rel¬ 
ativistic axions. The dominant contribution comes from 
the smallest loops at a given time which is some frac¬ 
tion a of the horizon f, a ~ k, where k — 0.15 is the 
backreaction scale (in units of t). The axion spectrum is 
dominated by the shortest wavelength ~ at. After the 
axion mass switches on at t\, strings transform into do¬ 
main walls at the moment t w « 1.55 1 1 , when further 
axion emission is terminated. Free axions became non- 
relativistic, with a conserved occupation number in each 
mode. Beginning at the moment t w , the net free axion 
field <f)(t, x) obeys the Sine-Gordon equation. As initial 
conditions at t = t w , we take <p(t,x) to be a superposi¬ 
tion of axion waves with uncorrelated phases, emitted by 
a number of loops, with the spectrum peaked at the mode 
k w ~ (atw)^ 1 . In the first approximation, <j>(t,x) can be 
modeled by a random Gaussian field with this spectrum. 
Since we will be interested in the parametric resonance 
at t > tQCD of the modes k which are many orders of 
magnitude larger than k ul , we will model the background 
field cj>(t, x) as a quasi-homogeneous field, with the scale 
of the inhomogeneity L ~ at w a(tQc d)/ a(t w ). The peaks 
of the random field 9o(x) can have a field value signifi¬ 
cantly exceeding the estimation (|23|). Since k w <C m a , for 
the epoch t > tQCD we can neglect the gradient term in 
the Sine-Gordon equation for the background field. At 
the scales L over which the net axion field <j>(t, x) can 
be considered to be homogeneous, it oscillates with the 
frequency m a . This is an approximation to the actual 
oscillations, which may be neither perfectly coherent nor 
perfectly harmonic. 

IV. AXION PARAMETRIC RESONANCE 
A. Axion Fluctuations 

We represent small fluctuations of the axion field by 
the modes (f>k{t)e lk x , where k is the comoving momen¬ 
tum. These fluctuations, (pkifij- obey the equation 

4>k + 3 Hcj) k + + to 2 cos e) <f> k = 0 . (24) 
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One might expect the parametric instability of the fluctu¬ 
ations due to the quasi-periodic variation of the effective 
frequency in Eq. (|24|). However, one should also wonder 
whether the expansion of the universe might suppress 
this effect pUl^l. 

The period of axion oscillations 

^ ~ 8 x 10” 10 (^ 51Q , k cV )s; meanwhile, the typical 

time for expansion of the Universe is — tQCD ~ 

7.7 x 10~ 6 s. Therefore, in the first approximation one 
can neglect the expansion of the universe in Eq. (|i|), and 
calculate the characteristic exponent ///,. for the growth 
of fluctuations, <pk oc e llkrnat . Then, one can take into 
account the adiabatic time dependence of and cal¬ 
culate the net effect as <f>k oc exp f t D dt fik(t) m a - We 
will see that this results in a good approximation to the 
actual evolution as determined by direct numerical solu¬ 
tion of the mode equation (|24|). 


B. Resonant Solution with Expansion 

We can now incorporate the expansion of the universe. 
Fortunately, this also can be done analytically, because 
the period of axion fluctuations ~ m” 1 is much smaller 
than the Hubble time H^ 1 around tQCD- We incorporate 
expansion into our static universe results by including the 
redshifting of the wavelength, k —> k/a , and the dilution 
of the background field, 9q —> d^/a 3 . Now, instead of ex¬ 
pression Q, we have the time-dependent characteristic 
exponent 


Adc(*) = 



4 k 2 

2 9 


(25) 


From here it follows that the expansion of the universe 
can shift the resonant modes out of the resonance. How¬ 
ever, as the restructuring of the band via expansion is 
somewhat faster than the redshifting of modes, the modes 
actually will be blueshifted relatively to the instability 
zone. That is, the upper (larger momentum) edge of the 
instability band will shift toward lower momentum val¬ 
ues faster than the modes are redshifted down. From ( |25| ) 
one can immediately estimate the time tk during which 


a given mode will be amplified, tk = ( e °^ a ) tQCD ■ In 
Fig. [| we plot the evolution of fluctuations as a function 
of time, obtained by a numerical integration of Eq. (24), 
which confirms the analytical estimation of tk- 



FIG. 4. The log of the axion fluctuations \(/)k\ 2 in 
an expanding universe as a function of time (in units of 
tQCD )• We plot the maximally excited mode, k « jr6 0 m a for 
9 0 = 0.1 (heavy line) in an expanding universe and the mode 
k = =9 0 m a that would be maximally excited in a static uni¬ 

verse (light line). As discussed in the text, the upper edge of 
the instability band shifts beneath these modes after a time 
tk = ( 9 ° 2 ”“ ) 4 tQCD, which is 42 x tQCD and 4 x tQCD for the 
modes shown. 


Amplification of fluctuations by the resonance is 
determined by the number of e-foldings: In 1| — 


ff k dt^k{t)m a , where we have assumed no back reac- 


JtQCD 

tion for the moment. From Eq. (25) we find 

rtk 


dt Hk(t) m a 


'tQCD 

= J%_ 


(Va k - 1 - arctanVafc - l) ( m a tQCD ), (26) 


where, ak is the change in scale factor from t-QCD to tk- 
Specifically, 


ak = 


I t k _ / Q 0 m o 
tQCD \ 2 k 


(27) 


Expression (]2(3|) reaches its maximal value 
ItQ CD dt MfcW ~ 0-091 (m a t QC D)d o for k « \0 o m a . 
Thus, the leading amplification of the amplitude of a fluc¬ 
tuation is given by 


<^ fc ~ ^Q0.039(m a tQCD)&Q 


(28) 


This formula is the main technical result of the paper 
regarding the axionic resonance. The strength of the res¬ 
onance depends on two numerical factors: the effective 
number of axion oscillations within a Hubble time around 
the QCD epoch, m a tQCD, and the value of the axion mis¬ 
alignment angle at tQCD, do- In Fig. || we plot the time 
evolution of the axion fluctuations obtained with the nu¬ 
merical solution of Eq. (J24| ) together with that derived 
with analytic formula (|26|) for jf QCD dtHk(t)m a for arbi¬ 
trary t. This demonstrates the reliability of our analytic 
anzats. 
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FIG. 5. The log of the axion fluctuations \(frk \ 2 in an ex¬ 
panding universe as a function of time (in units of tQCD )■ We 
plot the maximally excited mode, k ~ - 8 0 m a , for 0 o = 0.1 
(upper curves) and 0.05 (lower curves). The solid lines are 
from an exact numerical integration of the mode equation 
© in an expanding universe. The dotted curves are the 
approximation: ln|0fc| 2 « 2 J* dt fj,k(t)m a ', the integral is 
given by formula (|2^) evaluated for an arbitrary t > tQCD ■ 


C. Efficiency of the Axion Resonance 


Let us estimate the efficiency of the axion resonance in 
Eq. (|2^) for a viable range of parameters. For the axion 
mass range m a ~ 10~~ 3 — 10 -4 eV and tQCD — 7.7 x 10 -6 

sec, we have m a tQ C D ~ 5.85 x 10 6 ( 5 . 10 T 4 ° e y ) • Despite 


this large value for the number of axion oscillations in 
a Hubble time, the amplification of axion fluctuations 
is not generally very large because of the small average 
misalignment angle, 9q, estimated by formula (23). How¬ 
ever, as was mentioned at the beginning of this section, 
in most interesting cases the background field has large 
scale inhomogeneities. Therefore, we expect to have re¬ 
gions where 9q exceeds the estimation (|23|). In this case, 
only the r.m.s. value of 9q is given by (|23[). 

To amplify quantum axion fluctuations to the level 
where fluctuations are comparable with the background 
field (Sc/) 2 ) ~ 0o, the exponent in ( p8| ) should be as large 
as 0.039 ( m a tQCD ) 0 q > 23. Thus, regions of the universe 
where the background phase exceeds the level 


0q > 10 


-2 


Gr 


m a 


10- 4 eV 


- 1/2 


(29) 


will have significantly amplified quantum fluctuations of 
axions. Additionally, it has been noticed |l^| that non¬ 
linear axion dynamics prior to tQCD will further enhance 
the amplitude in the initially high 9 0 regions. If there 
are small scale, classical axion fluctuations, the thresh¬ 
old value of 9q for the effective resonance will be even 
lower. The wavelength of the primary resonant mode, 
k « ^9om a , is significantly smaller that the scale, L, 
over which 9q is expected to be inhomogeneous 


Lk ~ 10 2 9 0 t Q cDm a » 1 . (30) 

The region of excessive 9 0 corresponds to a proto¬ 
minicluster. The larger 9q is, the stronger the resonance 
effect will be. This means that in regions of particularly 
large 9 o, the background field will decay into axion fluc¬ 
tuations within several Hubble times ~ O(10)tQCD • The 
resulting axion fluctuations can be viewed as a super¬ 
position of very slowly moving (non-relativistic) axion 
waves of wavelength fc -1 . Thus, parametric resonance 
will cause the (quasi)homogeneous held in the region of 
a proto-minicluster to become strongly inhomogeneous. 
Then, when axions come to gravitationally dominant, a 
large number (~ (Lk) 3 ) of dense axion clumps will form 
from the remnants of the proto-minicluster. 

We have assumed that the background fluctuations in 
the region of proto-minicluster are harmonic and coher¬ 
ent. The actual fluctuations may not be perfectly co¬ 
herent. In this case, one would expect the resonance to 
be significant if 9q exceeds a level which is somewhat 
greater than the idealized estimation (p9|). To study the 
resonance in the case when the background fluctuations 
are not perfectly coherent, one could use lattice simula¬ 
tions like those in but with much greater resolution. 
Alternatively, one could try to mimic the background 
fluctuations as a random one dimensional process 0(t). 
These are, however, beyond the scope of the current pa¬ 
per. Also, we cannot extend our results to the scenarious 
where axions at tQCD are relativistic [fnj . For now, we 
consider our results as an indication of possible strong 
axion resonance in the regions of proto-miniclusters. 


V. PREHEATING IN NATURAL INFLATION 

Another important application of the theory Sine- 
Gordon parametric resonance elaborated above is to par¬ 
ticle production in the Natural Inflation scenario [||, 
where the pseudo-Nambu-Goldstone boson 0(f) plays the 
role of the inflaton field with the potential (ill). The 
typical parameters of Natural Inflation are / ~ M p . 
rri 0 ~ 10 12 GeV, and A ~ 10 15 GeV. First, we consider 
the rate of creation of inflaton quanta due to the self¬ 
interaction. This problem is reduced to the Sine-Gordon 
parametric resonance in an expanding universe. Then we 
consider other channels of inflaton decay. Although that 
issue goes beyond the stability of fluctuations in the Sine- 
Gordon model, we add this to the paper to complete the 
picture of reheating and preheating in the Natural Infla¬ 
tion scenario. 


A. Self-interacting Inflaton Decay 

The homogeneous inflaton field (j>(t) in this model 
obeys the equation 

0 + 3 Hcj) + m|/ sin j = 0 . 


(31) 








During the inflationary stage, the field is slowly 
rolling down to the minimum of the potential (jl|) . After 
inflation, the field oscillates about the minimum of the 
potential with a frequency m After several cycles, these 
oscillations can be simply approximated by harmonic 
ones, <j)(t) ~ <f>(f) sinm^t, with decreasing amplitude, 
<!>(£) = ^ 2 . The initial amplitude <f>o — O(l0~ 2 )M p . 
Inflaton oscillations correspond to the matter dominated 
equation of state for which a(f) a t 2 / 3 . 

Let us now turn to the vacuum inflaton fluctuations, 
5(f>. The temporal part of an eigenmode, </>*,(£), obeys the 
equation 

4> k + 3H<j) k + (^ + ™?<t> cos 4> k = 0 , (32) 


where the oscillating background phase 9 = ^jr- (c.f. 
Eq. (|24|)). At the start of oscillations, the typical ex¬ 
pansion time, H -1 , is comparable with the period of os¬ 
cillations and exceeds it at the latest stages. We find 
that the parametric resonance of inflaton fluctuations is 
insignificant. This result can be understood as follows. 
Resonance occurs in a single resonance band which spans 
momenta k/a(t) between zero and ~ 9o(t)m < i > , where 
6 0 {t) = = 0oa -3 / 2 . For this reason, resonant modes 

will not be redshifted out of the band, as was suggested 
in jj], Just the opposite occurs; momenta redshift more 
slowly than the background field amplitude, and they be¬ 
come blueshifted relative to the res onanc e band, just as 
occurred in the axion case, see Sec. IV B . Let us assume 
for a moment that the expansion is slow. Then we can 
use the approximation ( |25| ) for the characteristic expo¬ 
nent. The the number of e-foldings J dt will 

again be proportional to where to is the start 

of oscillations. However, the factor m^o is of order one, 
not 10 6 as it was in the axion case. Also, the phase 9 is 
necessarily rather small. The net result is, at best, a few 
e-foldings. Now, as the actual expansion at the begin¬ 
ning of the oscillatory phase is not slow, the resonance 
will be even less efficient. Thus, we confirm the claim of 
[|9i] that the parametric resonance of the inflaton field due 
to the self-interaction in Natural Inflation is insignificant, 
although our explanation is somewhat different. 


B. Resonant Inflation Decay with a g 2 <fr 2 X 2 
Interaction 

However, to make the picture of reheating in the Nat¬ 
ural Inflation scenario more comprehensive, we have to 
consider other channels of inflaton decay. Let us consider 
the decay of the inflaton into other bosons Xi via a cou¬ 
pling L g 2 4> 2 X 2 ■ Using a decomposition of the Bose field 
X into eigenfunctions Xk{t) e~ lk ' x with comoving momen¬ 
tum k, we obtain, for their temporal parts, 

Xk + 3 H\k + ^-^2 + = 0 ■ (33) 


where the back grou nd oscillations </>(£) are given by the 
solution of Eq. (|h]) . 

Eq. ([hi]) is an oscillator-like equation with a periodic 
effective frequency. Therefore, we should expect resonant 
solutions Xk ~ exp f dt/j,k(t)m ( j ) . Indeed, when the am¬ 
plitude of the background phase is sufficiently low, the 
background oscillations correspond to those in a massive 
theory -m 2 ^ 2 . Therefore, Eq. ( |33| ) is similar to the well- 
studying case of preheating in i?7i 2 i?!> 2 -inflation [Q. We 

introduce the usual resonance parameter q = 9 % . If 
q 1, we will have the creation of y-particles in the 
broad resonance regime. However, due to the expan¬ 
sion of the universe there will be no distinct resonance 
bands, but rather a broad range of resonant momenta 
Afc ~ q 1 / i m < i > where the resonance occurs in a stochas¬ 
tic manner pi. The growth of the particle’s occupation 
number rife (£) will be a stochastic process. For moderate 
values of q we approach the regular resonance regime. 



FIG. 6. The typical stochastic resonant production of 
particles at the particular choice of comoving momentum 
kg = 0.1 m|, and q = 10 2 with expansion. Time is mea¬ 
sured in the number of inflaton oscillations since the end of 
inflation. 


Let us estimate the range of allowed q values. From 
the basic parameters of the Natural Inflation scenario, 
we have q = — g 2 10 10 . The coupling parameter g 2 

,L <t> 

is restricted by the condition that radiative corrections 
do not alter the potential (|I]), Q^ig A <f) A log ^ < A 4 . This 
gives us g < -y- ~ 5 • 10 -4 , and finally q < 10 3 . In Fig. 
|(i| we show the evolution of the number density nk{t) of 
X'-particles calculated from a numerical solution to the 
equation for fluctuations (33). Thus, for the resonance 
parameter 1 < q < 10 3 the parametric resonance is very 
effective, and leads to the preheating of bosons %. The 
leading channel for inflaton energy transfer will be that 
with the greatest q. 

Additionally, it was found very recently that there can 
be preheating of fermions ip due to an inflaton coupling 
hipcj)^ [^2|. The creation of fermions occurs in the regime 
of parametric excitation, which is very different from the 
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perturbative calculations used earlier |jj for Natural In¬ 
flation. In perturbative reheating, the reheat tempera¬ 
ture was found to be rather low, T r ~ glO 8 GeV [||, 
where g 2 is coupling between inflatons and other parti¬ 
cles. However, as we have seen, the preheating of bosons 
and fermions in Natural Inflation can be very efficient. 
Therefore, there it is a possible to have a larger value 
of the reheating temperature in this scenario than that 
calculated with the perturbation theory. 

VI. SUMMARY 

In this paper, we have discussed the general theory of 
parametric resonance of fluctuations around an oscillat¬ 
ing homogeneous scalar field with a Sine-Gordon poten¬ 
tial. We found a closed form analytic resonance solution 
which shows that, in such a theory, there is a single insta¬ 
bility band which starts at zero momentum. The char¬ 
acteristic exponent and width of the resonance band are 
strongly dependent on the dimensionless amplitude, or 
misalignment angle, do, of the background oscillations, 
decreasing as do decreases. Although the background os¬ 
cillations are well described by a massive theory in the 
regime of small do, the correct properties of the para¬ 
metric resonance in this regime cannot be obtained by 
making the small angle approximation in the equation 
for fluctuations, which would lead to a Mathieu equa¬ 
tion. This is due to the non-perturbative nature of the 
parametric resonance. 

As a particular application of the general theory of 
Sine-Gordon parametric resonance, we considered the 
stability of axion fluctuations in the cosmological model 
with cosmic axions. By the QCD epoch, axions are 
non-relativistic and can be considered as a coherently 
oscillating background field with a misalignment angle 
do and frequency m a ~ 10~ 4 eV. The expansion of the 
universe is important to the strength of the resonance, 
but is not itself deadly as the typical expansion time sig¬ 
nificantly exceeds the period of axion fluctuations, e.g. 
iriatQCD 1. If the amplitude of fluctuations is as 

small as do < 1CV 3 , as required of the r.m.s. value 
in an axion dominated cosmology, then the expansion 
of the universe overdominates for viable ranges of m a 
and axionic parametric resonance is inefficient. How¬ 
ever, there are axion models for which the background 
field d 0 (x) at the QCD epoch is not perfectly homoge¬ 
neous. If there were no resonance, the regions of large 
do would become axion miniclusters when axions come 
to gravitationally dominant. Assuming the background 
oscillations are perfectly coherent in the regions of proto¬ 
miniclusters, we found that the axion fluctuations are un¬ 
stable in the proto-minicluster regions where do > 10~ 2 
for m a = 5 x 10~ 4 eV. The wavelengths of the reso¬ 
nant modes are significantly smaller than the expected 
scale of proto-minicluster regions. Thus, in these regions 
background axions at rest will resonantly decay into mov¬ 


ing non-relativistic axions. In other words, a smoothly 
varying axion field do(x) will decay into a strongly inho¬ 
mogeneous axion field. As a result, at the epoch when 
axions become gravitationally dominant, instead of a sin¬ 
gle, massive axion minicluster we will have a minicluster 
composed of a large number of smaller and denser axion 
“mini-miniclusters”. This may be important for axion 
detection strategies. The effect of axion resonance in the 
regions of miniclusters should be checked with lattice sim¬ 
ulations or with a more sophisticated analytic treatment 
of the background oscillations which may not be perfectly 
coherent. Therefore, we consider our calculations as an 
indication that there may be strong resonance effects in 
axion miniclusters. 

Another application of the general theory of the Sine- 
Gordon parametric resonance is to particle production 
in the Natural Inflation scenario, where the inflaton field 
has a Sine-Gordon self-interaction. We considered the 
production of inflaton fluctuations and found that the 
resonant decay of inflatons is not effective. Beyond the 
Sine-Gordon self-interaction, we consider the inflaton de¬ 
cay in this scenario via a four-leg interaction with other 
bosons. The channels of decay for which the coupling pa¬ 
rameter is in the region 10~ 5 < g < 5 • 10 -4 will lead to 
a broad parametric resonance; smaller g will correspond 
to a narrow resonance. Together with the preheating of 
fermions, this opens the possibility of preheating in the 
Natural Inflation scenario. 
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APPENDIX A: BACKGROUND SOLUTION OF 
THE SINE-GORDON EQUATION 

Here we derive the general homogeneous solution of 
the Sine-Gordon equation (|2:) for an arbitrary initial am¬ 
plitude <f) o, not necessarily much smaller than /, with no 
expansion of the Universe taken into account. Without 
loss of generality, we can restrict ourselves to the finite 
motion \<f>(t)\ < nf, with the constant total energy E 
given by (||). Let us put y(t) = Then from the 

energy integral ([|) one finds 

y 2 = ^(l+y 2 )( E -(2m 2 J 2 -E)y 2 ) . (34) 

Let us now recall the equation for the elliptic cosine, C = 
cn(T;IC), as a function of some variable r where /C is its 
modulus: 
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C 2 = 1C 12 + ( K 2 - K' 2 )C 2 - JC 2 C A , K' 2 = 1 - K 2 . (35) 


Here, () stands for the r—derivative. It is easy to 
check that Eq. ( |34| ) can be reduced to an equation 
for the elliptic cosine by means of the substitutions 
y = . 2 m 2 ~^~ cn ( r ; T = rn ^ where a dimension¬ 

less parameter is e = E/m 2 f 2 and the modulus is 
/C = yj e/2. Therefore, for the background oscillations we 
have tan -^^cn{m a t\ 1C), which finally results in 

Eq. (§• The (dimensionless) period of the background 
fluctuations, T t , is equal to 


s(l-e) 


T t = V 8 


da. 

a f) 


e — 1 + cos 


7T/2 

= Vs / , M 

o W (1 + sin 2 d) ^1 + (1 — e) sin 2 


(36) 


which is equal to 4K 



APPENDIX B: ALGEBRAIC FORM OF THE 
LAME EQUATION 


Here we demonstrate, for the sake of completeness, 
that the equation for fluctuations 4>k(f) i n the Sine- 
G or don model (0) (without expansion) can be reduced 
to the Lame equation with a particular choice of the co¬ 
efficients. Indeed, the general algebraic form of the Lame 
equation is pijj l 


Y" + 

+ 


1 


1 


1 


+ 


1 


2 \ z — ei z — e 2 z — e 3 
M + B Y = 0, 


Y‘ 


4 (z - e 1 )(z - e 2 )(z - e 3 ) 


(37) 


where Y(z) is the solution of the Lame equation. The 
equation for fluctuations (0) has the form (|37]) with e\ = 
1, e 2 = 1 — e, e 3 = —1, A = —2, n = 1, B = A = —2ft 2 . 
Therefore, <j>k(t) is a solution of the Lame equation with 
these parameters. Notice that the Lame equation has 
already been met in the context of parametric resonance 
in the theory of preheating with a XcjA effective potential 
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APPENDIX C: DERIVATION OF THE 
CHARACTERISTIC EXPONENT 


Here we show how one can derive Eq. ( 0 ) for the char¬ 
acteristic exponent /ifc from the analytic solution (14). 
Eq. ([[|) describes both solutions, 4>i(z) and (fi 2 (z). The 
resonant solution <j>(z) consists of four monotonic parts 


within a single period of the background oscillation. It 
turns out that at different quarters of the period either 
(f>i(z) or 4> 2 {z) corresponds to the exponentially grow¬ 
ing solution. Indeed, the square of the resonant solution 
within the first quarter of a period, while z is decreasing 
from 1 to (1 — e), is 

dz fi | Cfc ) 

M(z)\ y/(l — z 2 )(z — 1 + e) J 

(38) 

where M{z) is given by Cfc is given by 0), and 4>l 
is the square of the resonant solution in the beginning of 
the period when z = 1 . 

Within the second quarter of the period, while z is 
increasing from (1 — e) to 1, one has 


0 2 (z) = </>o exp 


0 2 ( 2 ) = <£ 1/4 ex P 


[ —(l- 

J M(z) V 

LL-e 


Cfc 


V(l-* 2 )(*-l + e) 

(39) 


where <j> 1/4 is the value of <f> z after the first quarter of the 
period = <f>(z = 1 ). 

Then the value of <fi 2 after half of a period is 


4>l/2 = 0o ex P 



dz 

M{z)Y/(J~l 2 ){z-l + e) 


(40) 


where the integral is understood as its principal value. 
The resonant solution has the generic form 0(z(r)) = 
P(z(r))e /iT , where P(z) is a periodic function. Since P 
has a period equal to half of the period of the inflaton 
oscillation, Eq. ( |40| ) is sufficient to find /i: 

hT t 
2 

1—e 


Cfc 


dz 


M(z)J(l-z 2 )(z-l + e) 


> 0. (41) 


The integral in this equation can be reduced to /(ft 2 ) 
given by (||). 


dz 


l-e 


M(z)y/(l-z 2 )(z-l + e) 


= 2 


7r/ 2 

f _ dd _ 

0 ^(1 + sin 2 d) ^1 + (1 — e) sin 2 d'j 


sin 2 d 


1 + ^1 — e + 2 ft 2 ^ sin 2 d 


= /(ft 2 ). 


(42) 
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